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Abstract. We study the tidal effects of a Kerr black hole on a neutron star in black 
hole-neutron star binary systems using a semi-analytical approach which describes 
the neutron star as a deformable ellipsoid. Relativistic effects on the neutron star self- 
gravity arc taken into account by employing a scalar potential resulting from relativistic 
stellar structure equations. We calculate quasi-equilibrium sequences of black hole- 
neutron star binaries, and the critical orbital separation at which the star is disrupted 
by the black hole tidal field: the latter quantity is of particular interest because when 
it is greater than the radius of the innermost stable circular orbit, a short gamma-ray 
burst scenario may develop. 



PACS numbers: 04.40.Dg; 97.60.Lf 
1. Introduction 

During the past decades, theorists have been modelhng various kinds of double compact 
objects since (1) they are among the most promising gravitational wave sources to be 
detected by ground-based and space-based laser interferometers [1] and (2) they have 
also been invoked as possible engines of short gamma-ray bursts [2] (see also [3, 4]) in 
the case of black hole-neutron star (BH-NS) and neutron star-neutron star (NS-NS) 
mergers. The remnants of both kinds of mergers, in fact, may result in a black hole 
with negligible baryon contamination along its polar symmetry axis and surrounded by 
a hot massive accretion disk: before the disk gas is accreted to the black hole, intense 
neutrino fluxes are emitted which, through energy transfer, trigger a high-entropy gas 
outflow off the surface of the accretion disk ("neutrino wind"); at the same time, energy 
deposition by vv annihilation in the baryon-free funnel around the rotation axis, powers 
relativistically expanding e^7 jets which can give rise to gamma-ray bursts [2]. The fate 
of BH-NS binaries, in particular, depends on the relative values of risco-, the radius of 
the innermost stable circular orbit, and rude^ the orbital separation at which the tidal 
disruption of the star by the BH occurs: if rjsco < ^tide the star is disrupted and then 
swallowed by the BH, otherwise it is swallowed without disruption. This is therefore 
a crucial issue for the SGRB mechanism we have described: only if rjsco < fude the 
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merger may result in the black hole + hot massive accretion disk + baryon-free funnel 
SGRB scenario. 

A significant effort in studying BH-NS coalescing systems has been undertaken 
in order to model them as gravitational wave sources and to understand gamma-ray 
burst engines. Unfortunately, BH-NS as well as BH-BH binaries (as opposed to NS- 
NS binaries) have not been observed yet and therefore it is not currently possible to 
infer their properties from observational data. Thus, in modelling the behaviour of such 
binaries and in making predictions about them, one has to rely on binary evolution 
and population synthesis models. Several approaches addressing different issues have 
been employed to study coalescing binaries. Post-Newtonian (PN) studies of compact 
binaries ([5] and references therein), for example, typically approximate the constituents 
of the binary as point sources. PN expansions may not converge rapidly enough in the 
strong-field region and thus are indicated for the inspiral phase: this has lead to a strong 
effort in stitching together PN studies with other methods that are more indicated for 
addressing the merging and ringdown, e.g. [6]. Other approaches deal with finite size 
effects due to at least one of the binary constituents; many of them assume Newtonian 
gravity in some or all aspects of the calculation. Finite size-effects are connected to tidal 
interactions between the binary constituents: these interactions are present well before 
the final phases of the coalescence and have been studied in the literature using various 
approaches and addressing several issues. We refer the reader to [7] for an overview on 
the literature. 

Nowadays, the mainstream strategy in studying BH-NS binaries is to develop fully 
relativistic codes, as is being done over the years by several groups (see [8], [9] for recent 
work). However the high computational cost of simulations of compact binaries makes 
it desirable to have frameworks to study certain features and dynamical regimes which 
admit the use of approximations, thus allowing a drastic reduction of the computational 
resources required. Motivated by this reason and by the intent of investigating SGRB 
engines, in this paper we focus on BH-NS binaries and develop a semi-relativistic model 
in order to describe the behaviour of a NS undergoing tidal interactions with its BH 
companion. 

Our starting point is the model developed in [10] by Wiggins and Lai: 

• it is essentially the affine model approach of Carter, Luminet and Marck [11], [12], 
i.e. the NS is treated as a Newtonian extended object which responds to its self- 
gravity, to its internal pressure forces and to the relativistic BH tidal field under 
the constraint that its shape is always that of an ellipsoid (hence the name "affine 
model"); this constraint descends from the fact that ellipsoidal figures are solutions 
to the dynamical equations describing a Newtonian incompressible star responding 
to the linear terms of a BH tidal field (see [13] for an extensive explanation). 

• NSs on equatorial circular orbits around a BH are considered and the quasi- 
equilibrium sequences of a Newtonian polytropic NS| interacting with a Kerr BH 

I Wiggins and Lai also consider white dwarfs, but we are not interested in white dwarfs in this paper. 
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are determined from large separation until the star tidal disruption: the critical 
orbital radius at which the tidal disruption occurs is compared to the radius of the 
ISCO in the context of SGRB engines. 

A similar approach has been used in [14]. 

In this paper we improve their model in two directions by 

(i) including general relativistic effects in the response of the NS to the BH tidal field 

(ii) using a general, barotropic equation of state (EOS) for the NS. 

Obviously, both points arise from the desire of building a model which has as many 
"realistic" features as possible; however there is more. The NS disruption, which may 
lead to the SGRB mechanism previously described, occurs when the BH tidal force 
starts prevailing on the NS self-gravity; this condition may be approximately expressed 
as aMj^s/RNS — ^bhRns I'^ude^ where a is a dimensionless coefficient, Af/vs and i?7V5 
are the NS mass and radius, and Mbh is the BH mass, rude hence depends essentially on 
two parameters: the mass ratio q = Mbh/Mns and the NS compactness at equilibrium 
C = Mns/Rns- The main advantage of Wiggins and Lai's model as opposed to many 
fully relativistic formulations, is that it allows one to choose large values of q. On the 
other hand, the weak point in ellipsoidal models in general, is that Newtonian self-gravity 
is adopted, and this choice is inappropriate to describe very compact stars. 

In our approach we improve this point by approximating the effects of general 
relativistic gravity by means of an effective scalar gravitational potential ^tov- This 
potential (see Appendix A) stems from the Tolman-Oppenheimer-Volkoff (TOV) stellar 
structure equations in General Relativity, and has been adopted and tested in the 
context of stellar core collapse and post-bounce evolution simulations: 

• in [15] Rampp and Janka present the Vertex code for supernova simulations; 
in order to approximate relativistic gravity, this code makes use of the generalised 
potential in place of the usual Newtonian potential in all Newtonian hydrodynamics 
equations 

• in [16] Liebendorfer et al. perform a comparison between results of the Vertex and 
Agile-BoltzTran codes, the latter being a fully relativistic (ID) hydrodynamics 
code; it is shown that both codes produce qualitatively very similar results except 
for some small (but growing) quantitative differences occurring in the late post- 
bounce evolution 

• in order to achieve a better agreement than that reported in [16], different 
improvements of the aforementioned effective relativistic potential are explored 
and tested by Marek and collaborators in [17]; the Newtonian equations of 
hydrodynamics remain untouched. 

The way we apply the main idea of [15]- [17] in the present, completely different, context 
is the following: the Eulerian hydrodynamics equations — which govern the star in the 
affine model and which are Newtonian — are left formally unchanged, but the pressure 
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profile of tlie star at equilibrium appearing in them is built with the TOV equations 
and gravitational potential used is ^tov- 

Furthermore, we extend the affine model to include general barotropic equations of 
state (see § 2). 

In this paper we concentrate on calculating rude with our approach. As a test of 
the model we reproduce the fully relativistic results obtained in [9], where the tidal 
disruption limit of binaries containing a Schwarzschild BH and a polytropic relativistic 
NS is calculated. We subsequently compute rude for binaries composed of a Kerr black 
hole and a NS, modeled with more realistic equations of state. 

The paper is organized as follows. In Section 2 we describe the basic equations of 
the affine model suitably modified as explained above; in Section 3 we test our model and 
present the results of the numerical simulations; in Section 4 we draw the conclusions. 

2. The Model 

In this section we build our BH-NS model whose essential features are: 

• the NS moves in the Kerr BH tidal field along the BH timelike geodesies 

• it maintains an ellipsoidal shape; more precisely it is a Riemann-S type ellipsoid, 
i.e. its spin and vorticity are parallel and their ratio is constant (see [13]) 

• its equilibrium structure is determined using the equations of General Relativity, 
its dynamical behaviour is governed by Newtonian hydrodynamics improved by the 
use of an effective relativistic potential 

• its EOS can be any (tabulated or analytic) barotropic EOS. 

In the spirit of [10], the equations for the NS will be written in the principal frame, 
i.e. the frame associated to the principal axes of the stellar ellipsoid, which is set up 
so that "1" denotes the direction along the axis that tends towards the BH, "2" is the 
direction along the other axis that lies in the orbital plane and "3" is associated with 
the direction of the axis orthogonal to the orbital plane. 

2.1. Overview of the Assumptions 

We shall assume that the star moves on an equatorial circular geodesic around the 
BH and neglect tidal effects on the orbital motion. We are thus working in the tidal 
approximation, i.e. we assume that: 

• Mbh ^ Mns (where Mbh is the black hole mass and M^s the star mass) 

• the stellar deformations do not influence its orbital motion. 

We neglect the perturbation that the star induces on the BH. 

According to the affine model, the surfaces of constant density inside the star form 
self-similar ellipsoids and the velocity of a fluid element is a linear function of the Xj's 
(the coordinates in the principal frame) [18]. This allows one to reduce the infinite 
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degrees of freedom of the stellar fluid to five dynamic variables and to deal with a 
(finite) set of ordinary differential equations governing the dynamics of the star. The 
five fluid variables are the three principal axes of the ellipsoid (oj, i = 1,2,3) and two 
angles and A which are defined by the differential relations 

d(j) ^ dX ^ /-,\ 

where r is the NS proper time, Q is the ellipsoid angular velocity measured in the 
parallel-transported frame associated with the star centre of mass and A characterises 
the internal fluid motion by 

aia2C 



A 



af + al 



(2) 



where ( is the (uniform) vorticity along the z-axis in the frame corotating with 
the ellipsoid. Given these assumptions, the Lagrangian governing the internal ("I") 
dynamics of the star may be written as 

Ci = Ct + Cb, (3) 

where "T" stands for "tidal" and "B" for "body" . In the next two sections we will write 
down both terms explicitly. 



2.2. BH-NS Tidal Interaction 
Ct may be written as 

Ct = - I d^Xp{x)<!)tide{x) 



1,...,3, 



(4) 



where Cjj and lij are the components of the BH tidal tensor and of the inertia tensor of 
the star. We recall that for a point mass orbiting a Kerr BH on an equatorial circular 
orbit (the NS centre of mass in our case) the energy and the ^-orbital angular momentum 
per unit mass are 

-2 - 2MBHr + a^/M^ 



E 



^/M^ir'^ - 2a^M^ + a 



where a is the black hole spin aligned with the orbital angular momentum, 

P = - ?>MBHr + 2a^/MBHr 
and out of E and one builds the constant 

K = {aE- L^f. 
The components of the tidal tensor for a Kerr spacetime are then 
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where r is the orbital separation, \1/ is an angle governed by 
E{L, -aE) + a 



(9) 



dr r^ + K 

which identifies the parallel-transporting frame associated with the star centre of mass 
[12] and (j) is the angle given in (1), that brings this frame into the principal frame by 
a rotation around an axis orthogonal to the orbital plane and passing through the star 
centre. 

The components of the tensor of inertia I appearing in (4) are defined as 

/« = /Apx.x,. (10) 

where p is the mass density distribution. In the principal frame, the tensor I is diagonal 
and takes the form 

I = Xldiagf^y, (11) 

where R^s is the isolated NS radius at (spherical) equilibrium, the a^'s indicate the 
lengths of the principal axes of the stellar ellipsoid and Ai is the star scalar quadrupole 
moment at spherical equilibrium (in isolation), i.e. 

M = -J x'x'pd^x = Y J ^^P^^ ■ 

Hereafter, the carets ( " ) denote variables referring to the isolated star at equilibrium. 
Notice that at equilibrium, as expected, the tensor of inertia and the scalar quadrupole 
moment are related by Tr(I) = 3M. 

2.3. Fluid Terms 

Following [10] we write 

Cb = Ti-U-V, (13) 

where T/ is the kinetic energy of the star internal (i.e. non-orbital) motion, U is the 
internal energy of the stellar fiuid and V is the star self-gravity potential. 
In general the internal kinetic energy of a body is given by 

Ti = j d'x^pu^, (14) 

where u is the velocity field of its internal motions which, for a Riemann-S type ellipsoid, 
is [13] u = Us + Ue, where 

Us = (^—A - X2ei + (^-—A + xica (15) 

is the spin velocity (the speed of the fiuid due to rotation), and 
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is the expansion/contraction velocity of the elhpsoid. Substituting these expressions in 
(14) we find 

2 



1 

+ 2 



(17) 



Since in this paper we are deahng with general barotropic EOS, the integrals appearing 
in Tj will have to be calculated numerically. 

The internal energy of the stellar fluid U is defined as the volume integral of the 
fluid energy density: 

U= I ed^x 



For any barotropic EOS, U is related to the pressure integral§ 



n = / P{p)d^x = 47r^^ P ( ^— ) r''dr (19) 

J ^NS Jo \aia2a3J 

by the following differential relation [11] 

dU = — dui] 
^ ai 

this equation is needed in order to write down the Lagrange equations. Finally, the 
self-gravity potential in general terms is defined as [13] 

du 



1 



(20) 



^ (af + n)A(n)' 

where 

A(m) = sj {al + u){al + M)(a| + m), 
V is the self-gravitational energy for the star at spherical equilibrium 

V = TrVij where V^j = - J dMdi{^)xj, (21) 

and $ is the gravitational potential. In spherical coordinates, the trace of Vij is 

i-Rns ^<|, 

V = -An / — r^pdf. (22) 

Jo df 

This expression must be evaluated numerically for the chosen barotropic EOS and 
gravitational potential $. In the case of a polytropic equation of state and for a 
Newtonian gravitational potential, the integral (22) can be performed analytically [13]. 
As already discussed, in this paper we improve the Newtonian modelling by using an 
effective gravitational potential in order to mimic relativistic effects in the framework 

§ Note that in the affinc approximation p = p/aia2a3, where p is the mass density for the star spherical 
equihbrium configuration. In this approach dM is unchanged while d^x — {aia2a3)/ R%gd^x. 
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of Newtonian hydrodynamics: equation (22) is where this potential must step in. In 
place of the Newtonian potential, following [15]- [17] who made the same assumption 
in hydro dynamical simulations of stellar core collapse and post-bounce evolution, we 
use the potential $tov which is solution of the equations of hydrostatic equilibrium in 
general relativity: 

d<^TOV _ [e(r) + P{r)][mTov{r) + Anr^Pjr)] 
dr p{r)r[r — 2mTov{i")] 

mTov{r)=ATT dr'r'^e{r'). (24) 
Jo 

Our goal is therefore achieved by plugging d^Tov/dr in (22) in place of d^/dr. Further 
details about this approximation are given in Appendix A. 

In Eqs. (23), (24) Mtov = f^Tov{R) is the gravitational (or TOV) mass of the NS, 
and the pressure (P), rest-mass density (p) and energy density (e) profiles are determined 
from the relativistic Tolman-Oppenheimer-Volkoff stellar structure equations. For the 
sake of simplicity, in these formulas we have omitted all carets on the equilibrium 
quantities. 



2.4- The Dynamics Equations 

The Lagrangian (3) for the star internal dynamics can easily be found by collecting the 
terms given in Eqs. (4), (11), (13), (17), (18) and (20); the Lagrange equations for the 
five fluid variables = {ai, 02, 03, 0, A} yield 

di = ai(A2 + _ 2a2An + -X^R^.a^A^ + - c^ai (25) 

27W M ai 

\ V ~ B? II 

02 = 02 (A^ + Q^) - 2aim + -—R\sa2A2 + 02202 

M ^2 

03 = - — i?;v5«3^3 + £3303 



M_ 

^NS 



Js — ^2 — Ci2{a\ — 



C =0 

where 11 is given by (19) and where we have defined 

du 

and 



A{u){a^ + u) 



Js = ^[{aj + al)n - 2aia2A], C = ^[{aj + al)A - 2aia2n]. (26) 

Js and C are, respectively, the spin angular momentum of the star, and a quantity 
proportional to its circulation in the locally nonrotating inertial frame. Notice that 
C = because we work in absence of viscosity. In this paper we consider irrotational 
models, i.e. with C = 0. 
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Equations (25) are a generalised version of Eqs. (31)-(35) in [10]; they are more 
general in the sense that (1) they are not restricted to the use of a polytropic EOS but 
are valid for any barotropic EOS and (2) they are written for any scalar gravitational 
potential for the self-gravity of the NS. When testing our programmes in order to 
reproduce the results of [10], we shall adopt the Newtonian potential in (22), i.e. 
d^/dr = 4:71 pr; elsewhere we will make the choice of using the effective TOV potential 



We now reduce the equations (25) to coupled algebraic equations by demanding 
quasi-equilibrium: during its evolution, in fact, a coalescing BH-NS binary will likely 
follow a quasi-equilibrium sequence with constant circulation (in particular, we set 
C = 0). Physically this statement relies on the fact that in such binaries the ratio 
between the orbital decay time due to gravity waves {tgw) and the tidal synchronization 
time (tsyn), i.e. the quantity governing the relative importance of viscosity, is smaller 
than unity [19]. This allows one to consider the fluid body not to be tidally locked 
with the orbital motion: internal fluid motion is therefore a necessary ingredient of the 
model. For a binary system in quasi-equilibrium one has to require that 



These are the equations we solve in this paper; in order to do so we adopt a 
Newton-Raphson scheme. The quasi-equilibrium sequence with constant circulation is 
parametrized by the binary orbital separation r. What we have is therefore essentially 
a "hydro without hydro" method, that is, in order to include finite-size effects in the 
inspiral phase we do not solve the hydrodynamic equations explicitly, but instead use 
"snapshots" generated by quasi-equilibrium conditions. To obtain a quasi-equilibrium 
sequence, we start by placing a non rotating spherical star in equilibrium at a distance 
To 3> Rns from the black hole||, we then gradually reduce the orbital separation, solve 
(27) and monitor the stars axes until a critical separation r^de is reached, at which no 
quasi-equilibrium configuration is possible. This critical distance is identified exploiting 
the fact that the algorithm cannot find any solution to the system (27) or by keeping 
track of the numerically calculated derivative drnorm/d{a2/ai), where 



which tends to zero at tidal disruption: both definitions yield the same values of rude^- 

II Wc make sure that the sequence we obtain is independent of tq. 

% An alternative approach to evaluate rude, based on the estimate of the Roche lobe radius in a 
Newtonian framework, has been used in [20]. 



(23)-(24). 



di = 
(f) = ^ 



(27) 





(28) 
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3. Results 

In this section we firstly compare our results on the tidal disruption of a TOV star to 
data extracted from the Full-GR results of Taniguchi et al. [9]: this allows us to assess 
the range of validity of our model and to prove the improvement gained over calculations 
performed with a Newtonian ellipsoidal star. Subsequently, we show and discuss some 
examples of quasi-equilibrium sequences in the case of a IAMq NS in tidal interaction 
with several kinds of stellar mass BHs; for the NS EOS we take two cases^: (1) the 
(stiff) EOS APR2, which yields Rns = 11.53km, and the (soft) EOS BALBNIHI, 
which instead gives Rns = 12.84 km. 

3.1. Comparison with Previous Results 

As a very first step we calculated BH-NS quasi-equilibrium sequences choosing a 
Newtonian gravitational potential for the NS, and a polytropic EOS with n = 1. We 
successfully reproduced the tables for rude given in Wiggins and Lai's work [10] and 
their data on equilibrium sequences. 

After these preliminary tests, we move on to using the improved gravitational 
potential for the NS and compare results given by our model with the Full-GR results 
presented in Figure 11 of [9] by Taniguchi, Baumgarte, Faber and Shapiro, who study 
the tidal disruption of = 1 and k, = 1 polytropic NSs, orbiting Schwarzschild BHs 
on equatorial circular orbits. In Figure 1 we show the results of this comparison; to 
facilitate the reader in analysing this figure along with Figure 11 of [9], we use the same 
physical quantities as Taniguchi et al. and hence the mass ratio 



is plotted function of 



which is the orbital angular velocity at which the tidal disruption occurs, normalised 
with respect to the polytropic length scale Rpoiy = y/^; in these definitions, Mbh is the 
hole mass (which is fixed by the value of q) and M^s is the ADM mass of the isolated 
star, which coincides with the gravitational-TOV mass. Of course, we take the same 
TOV equilibrium models used in [9], whose properties are discussed in the tables at the 
end of the aforementioned paper; the neutron star mass is fixed by choosing a specific 
value for the baryonic mass 



where rhTovi'f') is defined by (24). In Figure 1 we also show the results yielded by 
Wiggins and Lai's original approach: in this case the star model is fixed by equating the 

+ See 3.2 for details. 




(30) 




(31) 
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Newtonian star mass M^r to Mj^s- In each panel of the figure, the value of the stellar 
compactness C = M^s/Rns for the relativistic configurations is also displayed. 
The results displayed in Figure 1 show that: 

• our pseudo-relativist ic ellipsoidal model agrees with relativistic data much better 
than the Newtonian ellipsoidal model; 

• for a given compactness, the agreement between our data and Full-GR data 
improves as the mass ratio increases. This is due to the fact that our model assumes 
that the centre of mass of the star moves on a geodesic of the black hole spacetime 
in a reference frame centered on the BH, a condition which is better satisfied for 
larger mass ratios. For instance, for g = 9 and M^^ = 0.12Rpoiy {C = 0.1088), the 
results practically coincide; 

• for lower values of q the agreement between our data and the Full-GR data increases 
as the stellar compactness increases. 

For a more detailed and quantitative comparison, in Appendix B we tabulate the values 
of VttideRpoiy plotted in Figure 1. 
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Figure 1. Comparison between FuU-GR results from [9] (Full-GR), results obtained 
with the Newtonian ellipsoidal model of [10] (Newtonian) and results from our 
improved ellipsoidal model (Pseudo TOV). Each graph shows the mass ratio (Eq. (29)) 
versus the tidal disruption limit calculated as in Eq. (30). In all cases the neutron star 
EOS is an n = 1, K = 1 polytrope. The compactness C indicated in each panel refers 
to the Full-GR and to our models. 
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3.2. BH-NS Equilibrium Sequences and Tidal Disruption Limits 

Having tested the validity of our approach, we employ it to study different possible BH- 
NS binary configurations, determining the quasi-equilibrium sequences and the tidal 
disruption radius rude for equatorial circular orbits. We describe the fluid forming the 
NS with two different EOS proposed in recent years by the nuclear physics community, 
which we call APR2 and BALBNIHI. 

• The Akmal-Pandharipande-Ravenhall (APR2) hadronic EOS [21] describes matter 
consisting of neutrons, protons, electrons and muons in weak equilibrium; it is 
obtained within nuclear many-body theory using a variational approach to the 
Shrodinger equation; its microscopic input is based on the Argonne Wig potential 
for nucleon-nucleon interactions [22] — which is calibrated to deuteron properties 
and vacuum nucleon-nucleon phase shifts for laboratory energies Eiat up to 350 MeV 
— and on the Urbana IX [UIX) three-body potential [23]; relativistic corrections 
to both potentials are included [24] (which yields wig + 6v + UIX*) 

• The Balberg-Gal (BALBNIHI) EOS [25] describes matter consisting of neutrons, 
protons, electrons, muons and hyperons (S, A and H) in equilibrium. Assuming 
the mean field approximation, an effective potential is employed, whose parameters 
are tuned in order to reproduce the properties of nuclei and hypernuclei according 
to high energy experiments. This EOS is a generalization of the Lattimer-Swesty 
EOS, which does not include hyperons. 

APR2 is much stiffer than BALBNIHI and, in terms of stiffness, most modern EOS fall 
in between the two equations of state we choose. 

We consider several values of the mass ratio q = Mbh /M^Si ranging from 5 to 50, 
and three different values of the black hole spin: a = 0, a = Q.^Mbh and a = 0.99Mbh- 

In order to understand which sequences are possible candidates to be the engine 
of a SGRB, we compare the tidal disruption orbital separations rude with rjsco- In 
the Schwarzschild case, we determine the location of the ISCO through the analytic fit 
given in [9] , where the effects of the finite mass ratio and of the stellar compactness are 
taken into account: 

1 + q~^ 

rjsco = ^^^^6-3/2(l-0.444g-V4(i_3.54CV3))- (^2) 

Note that according to this equation the values of risco/ Mbh we give for a = in 
Figure 2 depend on the mass ratio q. 

Since no such fit exists for the Kerr BH-NS binaries, when a 7^ we estimate the 
ISCO by using the formulae derived in [26] for a point mass in the gravitational field of 
a Kerr BH: 

rjsco = Mbh{^ + ^2 T [(3 - Zi)(3 + Z^ + 2Z2)Y'^} (33) 
Zi = 1 + (1 - aVM|^)V3[(i + a/MsHf''' + (1 - a/MBHY'^] 
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where the upper sign holds for corotating orbits and the lower sign for counterotating 
orbits. Note that according to these equations the values that rjsco/^BH assumes for 
a 7^ do not depend on q. 

APR2 BALBNIHI 




. cl ^ 1 , 1 a_liJ , 1 , 1 ^ 1 , 1 n cl , 1 , 1 ^LJJ , 1 , 1 i 1 , 

3 4 5 6 7 8 9 10 3 4 5 6 7 8 9 10 




I gl L 1 , 1 , 1 , 1 , \ , 1 , 1 n 5I 1 ' 1 • 1 ' 1 ' 1 • 1 ' 1 

1 2 3 4 5 6 7 8 1 2 3 4 5 6 7 8 



Figure 2. 02/01 (continuous lines) and 03/01 (dotted lines) are plotted as functions 
of r/MsH for a NS orbiting a black hole. The angular momentum of the BH 
takes the values a/MsH = {0,0.5,0.99} as indicated in each panel. The mass ratio 
q = AIbh /Mpfs is indicated by the numbers next to each curve, which represent a 
quasi-equilibrium sequence. Mns = 1.4 in all graphs: in the left column the star 
is modeled using APR2 EOS (Rns = 12.84 km), in the right column using BALBNIHI 
EOS {Rns = 11.53 km). The ISCOs are shown by the dashed vertical lines, while the 
dot-dashed lines in the lower graphs mark the position of the black hole horizon. In 
the q = 50 case, the sequence does not terminate with tidal disruption; it stops when 
r = oi + R^fj^ where R%u is the size of the black hole outer horizon at the equatorial 
plane. 
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The results of our numerical integrations are displayed in Figure 2, where the 
ratios a2/ai (continuous lines) and fls/fli (dotted lines) among the NS axes are shown 
as functions of r /Mbh'-, the quasi-equilibrium sequences end when the tidal disruption 
of the NS is reached (see §2.4). The left panels refer to the EOS APR2, the right ones 
to BALBNIHI; the upper panels refer to the Schwarzschild case, the middle panels to 
a = 0.5Mbh, the lower panels to a = 0.99Mbh- The dashed, vertical lines are the 
locations of the ISCO. In the a = 0.99Mbh case we also indicate the location of the 
black hole outer horizon R'^j^ with a dot-dashed vertical line. From the graphs in Figure 
2 we extract the following information. 

• In each panel, the larger the mass ratio q is, the closer the NS can get to the 
BH before being disrupted, i.e. rude/MBH decreases as q increases. Therefore the 
conditions of the SGRB mechanism are more likely to be satisfied for low values of 

• If a = (first row) the star enters the ISCO before being disrupted for all the 
considered values of q, and for both the APR2 and the BALBNIHI EOS: therefore, 
SGRBs cannot be ignited. On the other hand, if a = 0.5Mbh (middle row) the 
star is disrupted at r > rjsco for q < 5.7 (APR2 EOS) and q < 6.6 (BALBNIHI 
EOS) and the SGRB may take place. If, finally, a = 0.99Mbh (third row), the star 
is disrupted at r > rjsco for g < 28 (APR2 EOS) and g < 33 (BALBNIHI EOS) 
and, again, the SGRB may take place. 

We note that as the the black hole spin increases, both rude/ Mbh and rjsco/MBH 
decrease; however, rjsco/MBH decreases more rapidly than rude/MBH, and 
consequently for a given mass ratio q, chances to develop an SGRB are higher 
if the black holes rotates faster. Thus, in general the conditions for the SGRB 
mechanism proposed in [2] to take place are favoured for low values of q and high 
values of a. 

• Comparing the two graphs in each row, we see that for NSs with a stiffer EOS, 
r tide/ Mbh is smaller. This is what one expects, since a star with a stiffer EOS 
is more compact and thus less deformable: hence a stronger gravitational field is 
needed to disrupt the star. 

• As the value of the BH spin increases, the values of the axes 02 (continuous lines) 
and 03 (dotted lines) tend to coincide. 

We also note that, if a = {].99Mbhi for g ^ 20 the star enters the ergosphere 
(r = 2M, since 6 = 7r/2) before disruption. 
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4. Conclusions 

In this paper we have studied the effects of a Kerr BH tidal field on a NS; we have focused 
on determining the orbital separation at which the NS is tidally disrupted {rude) and 
the star quasi-equilibrium sequence for several BH-NS binaries. Our work is based on 
the affine model, which we have improved with respect to previous works, by describing 
the NS self-gravity with an effective relativistic scalar potential (Eqs. (23) and (24)), 
and by using more realistic equations of state for the NS matter. This approach has 
the advantage of allowing a quick computation of rude and of the quasi-equilibrium 
sequences for any value of the binary mass ratio q and of the BH spin parameter a. 

A comparison of the results obtained using our approach for NS disruption with 
Full-GR results, which are available in the case of non-rotating BHs and polytropic 
NSs, shows an excellent agreement up to small values of q (which depend on the 
NS compactness), beyond which our model starts to underestimate rude] this happens 
because we assume that the centre of mass of the star moves on along a BH geodesic, 
a condition which is better satisfied for larger mass ratios. Our model agrees with 
relativistic calculations much better than ellipsoidal models in which the neutron star 
self-gravity is treated at a Newtonian level. 

Given the very good results of this test, we have determined the quasi-equilibrium 
sequences and the tidal disruption radii of several BH-NS binaries in order to evaluate 
with our model the role played by (1) the BH spin parameter a and (2) the equation of 
state of the NS fluid. We found that the soft gamma-ray burst scenario is favoured by 
low values of the mass ratio q and high values of the BH spin parameter a. In addition, 
the higher the BH spin parameter is, the more the values of the NS axes 02 and 03 (i.e. 
the principal axes which do not point at the BH) tend to coincide. 

In quantitative terms we have shown that, for the stiffer equation of state APR2 a 
SGRB may occur if g < 28 when a = 0.99 Mbh ot if q < 5.7 when a = 0.5 Mbh- For the 
softer equation of state BALBNIHI a SGRB may develop if g < 33 when a = 0.99 Mbh 
or if g < 6.6 when a = 0.5 Mbh- 

As a future development of our work, we intend to study the dynamical behaviour of 
BH-NS binaries and to determine the gravitational radiation emitted in the late phases 
of their inspiral. This means taking into account tidal interaction corrections to the 
orbital dynamics and to the gravitational waveform, and thus will require the inclusion 
in our model of a precise treatment of the orbit. Moreover we would like to improve the 
results of our model for low values of the mass ratio q. 

Appendix A. Effective General Relativistic Gravitational Potential 



In a Newtonian context, the gravitational potential of a non-rotating star at equilibrium 
is governed by the Poisson equation and the hydrostatic equilibrium equation, i.e. 

A<^Newt = 47rp (A.l) 



A Semi-relativistic Model for Tidal Interactions in BH-NS Coalescing Binaries 17 

d^Newt _ _ 1 dP 

dr p dr ^ 

and the system is closed by choosing an equation of state. 

In General Relativity, the equilibrium configuration of a spherical distribution of 
fluid is determined by the Tolman-Oppenheimer-Volkoff (TOV) equations of hydrostatic 
equilibrium 



dm{r) 
dr 



47rr^e(r) (A.2) 



dvir) ^m{r) + A-nr^Pir) 

dr r{r — 2m{r)) 

dP{r) (P(r) + e(r)) z/(r) 



dr 2 dr 

where m(r) is the gravitational mass enclosed within a radius r; also this system must 
be closed by choosing an equation of state. In order to define the effective relativistic or 
pseudo-relativistic scalar potential we use in this work, as in [15]-[17] we "mix together" 
equations (A.1)-(A.2), i.e. we describe the equilibrium configurations using the following 
equations 

dP (e + P) {rriTov + 4:7rr^P) 



dr r{r — 2mTov) 



(A.3) 



dr 

d^Tov 1 dP 

dr p dr 

where the first two are the first two TOV equations (A.2), and the last is the second 
Newtonian equation (A.l). The potential ^tov{'>^) obtained by integrating equations 
(A.3) is then used in the self-gravity tensor 



V,, = - J dM9,($)r, (A.4) 
which is needed to find V = Tr(Vij) used in the self-gravity potential (20). 

Appendix B. Tidal Disruption Data from Different Models 

In Table Bl, we provide the data plotted in Figure 1. A polytropic EOS with n = 1 and 
K = 1 is used to model the star. Each data set corresponds to one of the panels in Figure 
1, which is identified by the row indicating the NS compactness C and its gravitational 
mass normalised with respect to the polytropic length scale Rpoiy = k^^'^. The first 
column gives the mass ratio q. The remaining three columns are the orbital angular 
frequencies at tidal disruption, normalised with respect to Rpoiy/10, resulting from 
calculations performed with the three binary models considered: "Full GR" indicates 
data extracted from Figure 11 of [9], "Pseudo GR" data obtained with our model 
and "Newtonian" data obtained with the model of [10] in which the NS self-gravity 
is Newtonian. We remind the reader that the BH is always non rotating (the BHs in [9] 
actually have a small residual angular momentum, see the reference for further details). 
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Table Bl. Data plotted in Figure 1: each subtable corresponds to a panel of the 
figure, which is identified by the stellar compactness C and the gravitational mass 
normalised with respect to Rpoiy = k"/^. The mass ratio q is given in the first column, 
while the orbital angular frequency calculated at the tidal disruption limit resulting 
from each model considered here — Full-GR, Pseudo GR, Newtonian — is displayed 
in the form ^udeRpoiy / in the other three columns. See the text of the appendix for 
more details. 
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